A no-go theorem for a two-dimensional self-correcting quantum memory based on 

stabilizer codes 
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We study properties of stabilizer codes that permit a local description on a regular D-dimensional 
lattice. Specifically, we assume that the stabilizer group of a code (the gauge group for subsystem 
codes) can be generated by local Pauli operators such that the support of any generator is bounded 
by a hypercube of size O(l). Our first result concerns the optimal scaling of the distance d with 
the linear size of the lattice L. We prove an upper bound d = 0(L D_1 ) which is tight for D — 1, 2. 
This bound applies to both subspace and subsystem stabilizer codes. Secondly, we analyze the 
suitability of stabilizer codes for building a self-correcting quantum memory. Any stabilizer code with 
geometrically local generators can be naturally transformed to a local Hamiltonian penalizing states 
that violate the stabilizer condition. A degenerate ground-state of this Hamiltonian corresponds 
to the logical subspace of the code. We prove that for D — 1,2 the height of the energy barrier 
separating different logical states is upper bounded by a constant independent of the lattice size 
L. The same result holds if there are unused logical qubits that are treated as "gauge qubits". 
It demonstrates that a self-correcting quantum memory cannot be built using stabilizer codes in 
dimensions D = 1, 2. This result is in sharp contrast with the existence of a classical self-correcting 
memory in the form of a two-dimensional ferromagnet. Our results leave open the possibility for 
a self-correcting quantum memory based on 2D subsystem codes or on 3D subspace or subsystem 
codes. 

PACS numbers: 



I. INTRODUCTION 



One of the most intriguing open problems in the field 
of quantum information processing is whether one can 
build a self-correcting quantum memory — a macroscopic 
physical system storing quantum information for macro- 
scopic periods of time without active error correction [l[ . 
If it exists, such a hypothetical device could play the role 
of a "quantum hard drive" in a future quantum com- 
puting technology. It is also a question of fundamental 
interest how to design a system which allows for quantum 
coherence to be preserved at non-zero temperature, sim- 
ilar in some sense to the quest to understand the origin 
of high-temperature superconductivity. 

A possible scenario of how quantum self-correction can 
be achieved in realistic physical systems was proposed 
by Kitaev 0]. The main idea of [2j was to mimic self- 
correction in classical magnetic media storage devices us- 
ing quantum spin or electron Hamiltonians with topologi- 
cally ordered ground-states. Logical qubits encoded into 
the ground-state of such Hamiltonian become virtually 
isolated from the environment assuming that the envi- 
ronment can only probe the system locally by applying 
small static perturbations to the Hamiltonian. Note that 
what is usually called a topological order translates to 
coding language as simply the condition for a quantum 
code to have a macroscopic distance. A number of gener- 
alizations and proposals how this scenario could be imple- 
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mented in the lab have been made recently 0, 0, H, @, 0] • 
As was first explicitly noted in || , if the system inter- 
acts with a thermal bath, the presence of a topological 
order in the ground-state by itself does not guarantee 
self-correction. Indeed, one may expect that local errors 
created by thermal fluctuations will tend to accumulate 
which may result in a global error having a non-trivial ef- 
fect on the encoded qubits. This phenomenon of 'thermal 
fragility' of a topological order was studied by Nussinov 
and Ortiz in [§]. For Kitaev's 2D toric code the loss of 
topological order at any non-zero temperature was stud- 
ied m |a M, EE d2 . 

The main challenge in constructing a self-correcting 
memory is to combine the existence of a coding ground- 
space with a natural, passive, energy dissipation mech- 
anism which prevents errors from accumulating. It was 
argued by Dennis et al. @ that the 4D version of the 
toric code introduced by Kitaev in is indeed a self- 
correcting quantum memory. In [l[ Bacon considered 
the question of a self-correcting memory in dimensions 
less than 4. He introduced a three-dimensional model, 
based on a subsystem code which is now called the three- 
dimensional Bacon-Shor code, which could be a model of 
self-correcting quantum memory. 

In order to set the stage for our results on a possible 
self-correcting quantum memory, it is instructive to re- 
view the classical state of affairs. A classical ID Ising fer- 
romagnet has two degenerate ground-states correspond- 
ing to the classical repetition code. The distance of this 
code is n, the number of spins. A classical ID ferro- 
magnet is not a classical self-correcting memory, since 
the phase transition to a disordered phase occurs at zero 
temperature. The intuitive reason is that the energy cost 
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of a single domain of flipped spins is independent of the 
size of the domain and hence such domain can grow with- 
out cost destroying the order. A 2D classical Ising ferro- 
magnet is a self-correcting memory; the model exhibits 
a low-temperature memory phase separated by a phase 
transition to a disordered phase [HI- 

The ID and 2D Ising ferromagnets allow for storage 
of only classical information, since the distance of the 
classical repetition code with respect to phase-flip errors 
is 1. Guided by this example, one can take the following 
conditions as necessary for a classical spin system to be 
a self-correcting memory: 

1. The system has a degenerate ground-state; one has 
to flip a macroscopic number of spins in order to 
map one ground-state to another. 

2. A macroscopic energy barrier has to be traversed 
by any sequence of single-spin flips mapping one 
ground-state to another. 

Note that if condition 1 is violated, the environment can 
destroy the encoded information by acting only on a few 
spins. If condition 2 is violated, there is no reason to 
expect that an energy dissipation mechanism will prevent 
single-spin errors from accumulating into a logical error. 
For example, the ID Ising model satisfies condition 1 but 
fails to satisfy condition 2. On the other hand, the 2D 
Ising model satisfies both conditions since mapping one 
ground-state (all spins up) to the other (all spins down) 
requires creating a domain wall of macroscopic size. 

Let us now consider a system made up from n quantum 
spins (qubits) with a degenerate ground-state spanning a 
linear subspace C C (C 2 )® n which will be used to encode 
logical qubits [25| • The subspace C can be viewed as the 
logical subspace of a quantum code. 

In the present paper we shall restrict ourselves to sta- 
bilizer codes (a.k.a. additive codes). The main idea of 
stabilizer codes is to encode k logical qubits into n phys- 
ical qubits using a logical subspace C C (C 2 )®" spanned 
by states that are invariant under the action of a 
stabilizer group S, 

C = {\^j) g (C 2 )®" : P\i>) = |V>> VPg<S}. 

Here S may be an arbitrary Abelian subgroup of the 
Pauli group 

V = (iI,Xi, Zi,..., X n , Z n ) 

such that — / ^ S. For any stabilizer group S one can 
always choose a set of generators S = (Si,.. ., S m ) such 
that S a G V are self-adjoint Pauli operators. Thus the 
logical subspace C can be regarded as the ground-space 
of a Hamiltonian [26] 

m 

H = -^S a . (1) 

a=l 



Note that all terms in H pairwise commute. In addition, 
the Hamiltonian has a constant gap above the ground- 
space. All eigenvalues and eigenvectors of H can be ex- 
plicitly computed using the stabilizer formalism [l3|. If 
all generators are independent, one has m = n — k. In 
general m > n — k since it may be advantageous to choose 
an over-complete set of generators in order to make the 
energy barrier higher. 

We will assume that the physical qubits live at vertices 
of a D-dimensional lattice A = {1, . . . , L} D with open or 
periodic boundary conditions. This choice for regular 
D-dimensional lattice is made for pedagogical reasons, 
since all our results can be generalized straightforwardly 
to an arbitrary graph embedded into a ZJ-dimensional 
space. Accordingly, there are n = L D physical qubits. 
We would like the quantum memory Hamiltonian defined 
in Eq. (fT]) to be physically realizable, thus we demand 
that it involves only short-range interactions. More pre- 
cisely, we assume that the support of every generator S a 
can be bounded by a hypercube with r D vertices for some 
constant interaction range r. 

Examples of such Hamiltonians are the family of sur- 
face codes [3 , color codes [HI, , and 3D membrane- 
net models [17j . 

Recall that the weight \P\ of a Pauli operator P = 
Pi . . . P n G V is the number of non-identity single-qubit 
Pauli operators Pi . The distance d of a stabilizer code is 
the minimum weight of an undetectable Pauli error, i.e., 
an operator P G V preserving the logical subspace C and 
inducing a non-trivial transformation on C, 

d = min \P\. 

pec(s)\s 

Here C(S) is the centralizer of the subgroup S, i.e., a 
group of Pauli operators commuting with every element 
of S. The notation C(S)\S means a subset of elements 
in C(S) which are not in S. 

Let the system described by Hamiltonian Eq. fT} in- 
teract with a thermal bath. In order to decide whether 
the system is a good quantum memory, one in principle 
would have to choose a dynamical model describing the 
interaction with a thermal bath, then choose a particular 
decoding algorithm and calculate how the probability of 
failure at the decoding step scales with the system size. 
For realistic dynamical models such as Davies' weak cou- 
pling limit [HQS, 

pursuing this approach seems to be 
completely hopeless. However, it is not unreasonable to 
assume that the bath dephases the system on short time- 
scales T 2 <gc Ti Hence, after a short dephasing time, 
the state of the system is a mixture of energy eigenstates 
and the effect of the bath is a process of energy exchange 
with the system. The state of the system can thus be rep- 
resented as a probabilistic mixture of energy eigenstates 
E\ip), where G £ is the encoded logical state and 
E G V is a Pauli error. An error Be? destroys encoded 
information if it is undetectable by the code C, i.e., if the 
restriction of E onto C induces a non-trivial transforma- 
tion on C. Under these assumptions we can now state 
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the quantum analogue of conditions 1,2 which are nec- 
essary for a quantum spin system to be a self-correcting 
quantum memory: 

lq. The ground-space of a (Hamiltonian) system is the 
logical subspace of a quantum error correcting code 
with macroscopic distance. 

2q. A macroscopic energy barrier has to be traversed by 
any sequence of single-qubit Pauli errors resulting 
in an undetectable error. 

The question that we address in this paper is: what 
quantum error correcting codes are compatible with con- 
ditions lq and 2q? Note again that our conditions are 
necessary conditions for a self-correcting quantum mem- 
ory; we do not claim that they are also sufficient condi- 
tions. In particular, whether a model is a self-correcting 
quantum memory may depend on details of the scaling 
of entropy versus energy at non-zero temperature. 

A. Bounds on The Distance 

Our first result concerns the optimal scaling of the dis- 
tance d with the linear size of the lattice L. It allows one 
to check whether a code is compatible with condition lq. 
In Section [II] we prove the following theorem 

Theorem 1. Let S = (Si, . . . , S m ) be a stabilizer code 
on a D-dimensional lattice A = {1, . . . , L} D . Suppose 
the support of any generator S a can be bounded by a hy- 
percube with r D vertices. Then the distance of S satisfies 

d<rL D -\ (2) 

This bound holds for both periodic and open boundary 
conditions. 

This theorem implies that any stabilizer code on a ID 
lattice fails to satisfy condition lq. It contrasts the fact 
that for classical ID codes the distance can be of order 
L (consider as an example the ID repetition code). For 
a 2D lattice the bound Eq. (|2|) allows d to grow linearly 
with L which is compatible with condition lq. Surface 
codes [HI provide an example of 2D codes for which the 
bound Eq. ^ is saturated. 

We are not aware of any codes saturating the bound 
Eq. ^ for D > 3. Note that a generalized toric code on 
a Z?-dimensional lattice [28[ has distance d ~ lL-°/ 2 J . We 
conjecture that in D = 3 the distance of any stabilizer 
code with local generators satisfies d — O(L). 

If a stabilizer code has more than one logical qubit, 
some of the logical qubits may be protected from the en- 
vironment better than the others. Note that Theorem [1] 
bounds the minimum weight of logical operators for the 
worst choice of a logical qubit. One can ask whether 
the same bound applies to the best logical qubit as well. 
To state this question more formally assume that we en- 
code k qubits using a stabilizer code S with g + k logical 



qubits for some g > 0. We will regard the extra g logi- 
cal qubits with the corresponding logical Pauli operators 
X\, Zi, . . . , X g , Z g as unused "gauge qubits" [2(| such 
that any error affecting only the gauge qubits can be ig- 
nored. At the same time we can use the logical Pauli 
operators on the gauge qubits to minimize the weight of 
'useful' logical operators. Thus, in this case the relevant 
distance would be the minimum weight of a Pauli oper- 
ator P preserving the logical subspace, i.e., P S C(S), 
that cannot be generated by stabilizers and logical Pauli 
operators on the gauge qubits, that is, 

d(G)= min |P|, g=(S,X 1 ,Z 1 ,...,X g ,Zg). 
pec(s)\g 

(3) 

In Section IIII CI we will use the formalism of subsystem 
codes to prove the following extension of Theorem [1] 

Theorem 1*. Under the assumptions of Theorem]]] the 

distance d(Q) satisfies the bound Eq. 0) for any choice of 
logical operators Xi, Z\, . . . ,X g , Z g on the gauge qubits. 

Let us remark that an analogous bound does not apply 
to 'bare' logical operators, i.e., P G C(Q)\Q. The minimal 
weight of 'bare' logical operators can be of order L D , see 

B. Bounds on The Energy Barrier 

Our second result concerns the scaling of the energy 
barrier with lattice size L. In contrast to the distance d 
the energy barrier is not a property of the stabilizer group 
S only — it depends on the choice of generators S a used 
to define the Hamiltonian Eq. |T]). As was mentioned 
above, the set of generators can be vastly overcomplete 
(consider as an example the 4D toric code). In order 
to exclude overcomplete generating sets in which some 
generator S a appears a macroscopic number of times we 
shall impose a (natural) restriction that any qubit can 
participate only in a constant number of generators S a . 

Note that the Hamiltonian Eq. ([T|) has a ground-state 
energy —to. Let \ip) S £ be any ground-state and E G V 
be any Pauli operator. Since E either commutes or anti- 
commutes with every term in the Hamiltonian, a state 
E \tp) has energy (rp\E^ H E\ij)) = —m + e(E), where e(E) 
is proportional to the number of Pauli operators S a in 
the Hamiltonian H which anticommute with E, 

e(E) = 2#{a : S a E = -ES a }. (4) 

We shall refer to e(E) as the energy cost of a Pauli oper- 
ator E. 

Let us say that a sequence Pq, Pi, . . . , Pf G V is a walk 
on the Pauli group starting at Pq and arriving at Pt iff 
for all i the operators P, and Pj+i differ by at most one 
qubit. Let W(S,T) be a set of all walks starting at S 
and arriving at T. For any walk 7 G W(S, T) let t m ax{l) 
be the maximum energy reached by 7 

e-raaxil) = maxe(P). 

P£7 
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Suppose the environment tries to implement a Pauli error 
E G V by a sequence of single-qubit Pauli errors. The 
minimum amount of energy the environment has to inject 
into the system in order to implement E is 

d x (E) = min e max {l)- 

7GW(/,B) 

Thus the energy barrier the environment has to overcome 
in order to implement a non-trivial logical operator on 
the encoded qubits is 



d* = 



mm mm e m ax 

E£C{S)\S j£W(I,E) 



(7). 



(5) 



Let us emphasize once more that in contrast to the dis- 
tance d the energy barrier gF is not uniquely determined 
by the stabilizer group S but is a function of the gener- 
ating set Si,..., S, n used to construct the Hamiltonian 
Eq. Q. 

Note that our assumption that the environment im- 
plements a logical operator by single-qubit Pauli errors 
is not a restriction if we are interested in determining 
the scaling of the energy barrier gP with L. To see this, 
imagine that the environment instead makes Pauli errors 
with weight at most w = 0(1) at the time. One can 
always simulate a single step of such a walk by w steps 
with single-qubit Pauli errors. It can increase the maxi- 
mum energy reached by any walk by at most a constant 
and hence the energy barrier d x can change at most by a 
constant. 

Using the assumption that any qubit participates in 
0(1) generators S a one can easily prove a naive upper 
bound 



d t = 0{d). 



(6) 



Therefore, if we want cF to grow with L we must look for 
codes with the distance d growing with L. For 2D codes 
one may have d ~ L and thus the naive upper bound 
Eq. ([5]) cannot rule out a possibility that cP grows with 
L. In Section [IJ we prove the following theorem: 

Theorem 2. Let S = (Si, . . . , S m ) be a stabilizer code 
with local generators on a 2D lattice such that each qubit 
participates in a constant number of generators. Then 
the energy barrier eP is upper bounded by a constant in- 
dependent of the lattice size L. 

This theorem tells us that 2D stabilizer codes can- 
not be compatible with condition 2q and thus the cor- 
responding Hamiltonians defined in Eq. |T]) cannot be 
a self-correcting quantum memory. This result contrasts 
the fact that 2D classical codes with local generators (e.g. 
the 2D Ising model) can be used to build a self-correcting 
classical memory. We remark that Theorem [5] applies to 
both open and periodic boundary conditions. It can also 
be easily generalized to quasi-2D lattices with the num- 
ber of 2D layers bounded by a constant independent of 
L. 

Theorem [5] can be strengthened for stabilizer codes 
which encode more than one logical qubit. Specifically, 



we can encode k qubits using a stabilizer code with g + k 
logical qubits by regarding the extra g logical qubits as 
"gauge qubits" , see the discussion after Theorem Q] Ac- 
cordingly, any error affecting only the gauge qubits can be 
ignored. In this case the energy barrier that the environ- 
ment has to overcome in order to implement a non-trivial 
logical operator would be 

d t (Q) = min min e ma x{l), (7) 
Eec(s)\g T ew(z,fi) 

where Q = (S, X\, Z±,..., X g , Z g ) is a group generated 
by S and the logical Pauli operators on the gauge qubits. 
In Section IIII CI we will use the formalism of subsystem 
codes to prove the following extension of Theorem 

Theorem 2*. Under the assumptions of Theorem&the 
energy barrier d x (G) is upper bounded by a constant in- 
dependent of the lattice size L for any choice of the gauge 
Pauli operators X\,Z\,..., X g , Z g . 

It has been observed by many authors [l|, [l9[ that the 
2D toric code does not feature a macroscopic energy bar- 
rier since the logical operators for this code have a string- 
like geometry. The energy cost of a partially implemented 
logical operator comes only from the two end-points of a 
string which cannot be larger than some small constant. 
One can ask whether logical operators for general 2D sta- 
bilizer codes can always be chosen as "strings" (in which 
case the logical operators can be interpreted as moving 
point-like anyons around the lattice). The proof of The- 
orem [2] presented in Section |TT] provides a partial answer 
to this question. It implies that for any 2D stabilizer 
code with generators of size r x r there exists at least 
one logical operator whose support can be covered by a 
rectangle of size r x L. Such a rectangle can be regarded 
as a quasi-lD string. 



C. Beyond Subspace Stabilizer Codes 

Hamiltonians of the form in Eq. |T]) may seem like 
overly restrictive models to consider for a quantum mem- 
ory, in particular since they involve a set of commut- 
ing operators. Common models in many-body physics 
involve geometrically-local sets of non-commuting oper- 
ators, such as the Heisenberg model and its variants, 
Kitaev's honeycomb model [J], the quantum compass 
model, see e.g. [21] . etc. Such Hamiltonians can always 
be written as 



H 



E 



(8) 



where r a are real coefficients and G a are elements of the 
Pauli group V with local support. The formalism of sta- 
bilizer subsystem codes 20] provides a systematic way 
of constructing Hamiltonians of the form Eq. ([5]) with 
a degenerate ground-state that can be taken as possible 
models of a quantum memory. This formalism focuses 
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on those symmetries of the Hamiltonian that can be de- 
scribed using the Pauli group ignoring all other non-Pauli 
symmetries such as the S[/(2)-syrnmetry (since they can- 
not be analyzed using the framework of stabilizer codes). 

By analogy with stabilizer codes which can be fully 
characterized by the stabilizer group, a subsystem code 
can be characterized by its gauge group. In our case 
the gauge group Q is the group generated by the Pauli 
operators G a , see Eq. ([8]) . Conversely, any (non-abelian) 
subgroup Q C V with local generators Gi, . . . , G m yields 
a family of Hamiltonians Eq. {8} . The symmetries of the 
Hamiltonian H defined in Eq. © are described by the 
centralizcr of Q, i.e. those elements in V which commute 
with every element in Q. It can be always represented as 

C(G) = (S,X j ,Z j ,j = l,...,k), 

where S is the center of the gauge group, i.e. S = 
G nC(<5), and Xj,Zj are logical Pauli operators of the 
quantum code. Since the symmetries of H include the 
logical Pauli group on the k logical qubits, any eigen- 
value of H has degeneracy at least 2 fe . In particular, the 
ground-state of H can be used to encode k logical qubits. 

The group S is called a stabilizer group of the subsys- 
tem code. It induces a decomposition of the Hilbert space 
of n qubits into (^-invariant sectors, (C 2 )®" = © s £ s , 
where s labels different syndromes, i.e., irreps of the sta- 
bilizer group. Given any syndrome s, the corresponding 
sector C = C s defines a logical subspace of the subsystem 
code |3(J. The logical subspace possesses a subsystem 
structure, C = £i og icai ® C gauge , where Ci ogica i defines 
the logical subsystem and C gauge defines the gauge sub- 
system. The gauge group Q restricted to C can be iden- 
tified with the group of Pauli operators on the gauge 
subsystem C gauge . The group generated by the logical 
operators Xj , Zj restricted to C can be identified with 
the group of Pauli operators on the logical subsystem 
Ci og icai- Accordingly, dim (Ci ogica i) = 2 k . 

It is important to note that the stabilizer group of a 
subsystem code with local generators does not necessar- 
ily have local generators. This fact may make such sub- 
system codes potentially more powerful than stabilizer 
codes. Subsystem codes with an Abelian gauge group 
have no gauge qubits, i.e., C — Ci og i ca i hence they coin- 
cide with the (subspace) stabilizer codes discussed above. 
Note that for (subspace) stabilizer codes the stabilizer 
group S coincides with the gauge group Q. 

The distance d of a subsystem code is defined as the 
minimum weight of a Pauli error P G V preserving C and 
inducing a non-trivial transformation on Ciogicai, 

d= min [PL S = gnC(G). (9) 

PeC(5)\e 

Let us mention a useful identity 

C{S) = Q-C{G) (10) 

that holds for any subsystem code Q. Note that a subsys- 
tem code with a large distance must have a non-trivial 



stabilizer group. Indeed, if S = (I), one has two possi- 
bilities: (i) any single-qubit Pauli operator belongs to G, 
and (ii) some single-qubit Pauli operator is a non-trivial 
logical operator. In case (i) one has G — V, i.e., the code 
has no logical qubits. In case (ii) one has d = 1. 

The best known example of a locally generated sub- 
system code is the Bacon-Shor code (T], [22[ which can be 
defined on 2D or 3D lattice. The Hamiltonian H asso- 
ciated with the 2D Bacon-Shor code corresponds to the 
quantum compass model in condensed-matter physics, 

Hbs — — JxXijXi+i.j + J z ZijZij + i. 

i<i,j<L 

Here each pair i,j represents a vertex of the lattice and 
Jx, Jz are arbitrary real coefficients. Accordingly, the 
gauge group G of the 2D Bacon-Shor code has genera- 
tors XijXi+ij and ZijZij+\. This code has distance 
L. The stabilizer group S is generated by operators 
Sf = U.LiX id X i+lij and SI = II' : Z,,Z,„. : . One 
can easily check that all elements of S have weight at 
least 2L, that is, S cannot have local generators. The 
code has one logical qubit with logical Pauli operators 

A = n^iA M andZ = nii^,i- 

The example of the Bacon-Shor code demonstrates 
that some subsystem codes with local generators origi- 
nate from subspace codes with highly non-local stabilizer 
groups. It might suggest that subsystem codes can beat 
the upper bound Eq. @. We will show that this intuition 
is wrong by proving 

Theorem 3. Let G = (G\, . . . , G m ) be the gauge group 
of a subsystem stabilizer code on a D -dimensional lattice 
A = {1, . . . , L} D . Suppose the support of any generator 
G a can be bounded by a hypercube with r D vertices. Then 
the distance of G satisfies 

d<3rL D -\ (11) 

Note that the distance of a subsystem code depends 
only on the gauge group G- Therefore, condition lq is 
true or false for any Hamiltonian in the family Eq. ([8]) 
regardless of the choice of coefficients r a . 

As far as condition 2q is concerned, the first natural 
question is how to define the energy barrier dr that must 
be traversed in order to implement a logical operator. 
In Section IIII Bl we propose a definition of (ft that de- 
pends on the energy spectrum of the Hamiltonian in dif- 
ferent sectors C s . We also derive a usable upper bound 
on the energy barrier that does not require computing 
the energy spectrum. For any particular Hamiltonian 
this upper bound can be calculated using only the sta- 
bilizer formalism. Unfortunately, we have not yet been 
able to prove an analogue of Theorem [5] for subsystem 
codes. In Section llllBl we prove a very weak upper bound 
on (ft, namely, eP = 0(d). This leaves as an open ques- 
tion whether Hamiltonian models based on 2D subsystem 
codes may have a macroscopic energy barrier. 
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D. Organization of The Paper 

The proof of Theorems [I][3] relies on three technical re- 
sults, Lemma's [T][3] The first Lemma Q] which we call 
a "cleaning lemma" asserts that for any logical opera- 
tor P G C(S) and for any subset of qubits M, \M\ < d, 
one can choose a stabilizer S G S such that PS acts 
trivially on M. The stabilizer S involves only those gen- 
erators S a whose support overlaps with M. Lemma [5] 
proves an analogue of this cleaning lemma for subsystem 
codes. The third Lemma [3] which we call a "restriction 
lemma" relates the distance of a code S defined on the 
entire lattice A and the distance of a code Sm obtained 
by restricting generators of S onto some subset of qubits 
M G A. The restriction lemma applies to both subspace 
and subsystem stabilizer codes. 

Stabilizer codes with local generators are discussed in 
Section [IT] which proves Theorems [T] and [2] Subsystem 
codes with local generators are discussed in Section Hill 
We prove the upper bound on the distance of subsystem 
codes (Theorem [3]) in Section MI Al Upper bounds on 
the energy barrier of subsystem codes are discussed in 
Section fill B I Finally Section fill CI shows how to prove 
Theorems 1* and 2*. Some open problems are discussed 
in Section HVl 

II. STABILIZER CODES 

Recall that V denotes the Pauli group on n qubits. We 
assume that the qubits live at vertices of a D-dimensional 
lattice A = {1, . . . , L} D . For any subgroup S C V and 
any subset MCA introduce a group 

S(M) = {P eS : Supp(P) C M} 

which includes all elements of S whose support is con- 
tained in M. In particular, V(M) is a group of all Pauli 
operators whose support is contained in M. Introduce 
also a group 

S M = {Pe V{M) : PQeS for some Q G V(A\M)} 

which includes all Pauli operators P G V(M) that can be 
extended to some element of S. In other words Sm is a 
group obtained by restricting elements in S to M. Note 
that if some element P G S crosses the boundary of M 
then the restriction of P onto M is no longer element of 
S. By definition S(M) C S M C V[M) C V. 

Our main technical tool will be the following "clean- 
ing" lemma. It allows one to clean out any region M C A 
of size smaller than the distance such that no logical 
operator of the code contains X, Y, or Z on qubits of 
M. More formally, one can multiply any logical operator 
P G C(S) by a stabilizer S G S such that PS acts triv- 
ially on M. The stabilizer S uses only those generators S a 
whose support overlaps with M, see Figure [TJ The clean- 
ing lemma is particularly useful when the generators of S 
are local. In this case the cleaning changes P only inside 




FIG. 1: Support of a logical operator P G C(S) is shown 
by the colored dots (representing X, Y, and Z operators). In 
order to clean out the region M (interior of the dashed ellipse) 
one can multiply P with a stabilizer S G S. The stabilizer 
S includes only those generators S a whose support overlaps 
with M. It yields an equivalent logical operator PS acting 
trivially on M. The cleaning can performed simultaneously 
for different logical operators and for different spacial regions. 

M and in a small neighborhood of the boundary of M. 
Thus the cleaning of P can be done multiple times, so 
that multiple 'holes' can be made into the support of P. 

Lemma 1 (Cleaning Lemma). Let S — (Si, . . . , S m ) be 
a stabilizer code and M C A be an arbitrary subset of 
qubits. Denote J(M) a set of indexes a such that the 
support of S a overlaps with M . Then one of the following 
is true: 

(1) There exists a non-trivial logical operator P G C(S)\S 
whose support is contained in M , 

(2) For any logical operator P G C(S) one can choose a 
stabilizer 

S= II s a a > ^e{0,i} 

aeJ(M) 

such that PS acts trivially on qubits of M . 

Proof. Let (il) C V be a subgroup of phase factors. By 
definition of the subgroups Sm and S{M) one has the 
following inclusion: 

(il) ■ S(M) C C(S M ) n V(M). (12) 

If the inclusion Eq. (fT2")) is strict, there exists a Pauli 
operator P G V(M) such that P G C{S M ) but P £ S, 
that is, P G C(S)\S. It corresponds to the case (1) in 
the statement of the lemma. Otherwise, Eq. (fT2")) is an 
equality, that is, C(S M ) H V{M) = (il) ■ «S(M). Taking 
the centralizer of both sides one gets 

S M = C(S(M))r\V(M). (13) 

(Note that for any subgroup T of the Pauli group one has 
C(C(T)) = (il)-T.) Let P G C(S) be any logical operator 
and Pm be a restriction of P onto M. By definition, 
P M £ C{S{M)) n V(M). It follows from Eq. (T5§ that 
Pm G <Smj that is, one can find a stabilizer S G S such 
that PS acts trivially on M. Let us expand the stabilizer 
S in term of the generators and take out all generators 
whose support does not overlap with M. It yields a new 
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stabilizer S' G S such that PS' still acts trivially on M. 
On the other hand, S' is a product of generators S a , 
a G J(M) 1 so we arrive at the case (2) in the statement 
of the lemma. □ 

The proof of both Theorems [1] and [2] goes by applying 
the cleaning lemma to get an upper bound on a linear 
distance of a code that we define below. Let us first 
define a linear distance of a code when the lattice A = 
{1, . . . , L} D has open boundary conditions. 

Definition 1. Given a Pauli operator P G V define 
its linear dimension d\(P) as the minimum length of 
a contiguous interval R C {1, ...,L} such that R x 
{1, . . . , L} D ~ 1 covers the support of P. Given a stabi- 
lizer code S define a linear distance of S as 

d 1 (S)= mm di(P). (14) 
Pec(s)\s 

Thus a linear distance of a code characterizes the 
minimal linear dimension of non-trivial logical operators 
along some fixed coordinate axis. Above we have chosen 
the first coordinate axis but one could similarly define 
a linear distance for any other coordinate axis. For pe- 
riodic boundary conditions we define a linear dimension 
d±(P) as the minimum length of a contiguous interval 
R C Zl such that R x (Zl) d_1 covers the support of P. 
Then a linear distance of a code is defined by the same 
formula Eq. (TH)) . 

Proposition 1. LetS = (Si, . . . , S m ) be a stabilizer code 
on a lattice A — {1, . . . , L} D with open or periodic bound- 
ary conditions. Suppose the support of any generator S a 
can be covered by a hypercube with r D vertices. Then 

di{S)<r (15) 

for any L > 2(r — l) 2 . 

Proof. We shall prove the proposition for D = 2. The 
generalization to other dimensions is straightforward. 
Let us assume that the bound Eq. (fT5|) is not satisfied, 
that is, 

di(<S)>r (16) 

and show that it leads to a contradiction. The idea behind 
the proof is depicted in Fig. [5] 

One can easily check that any integer L > 2(r — l) 2 can 
be represented as L = a(r—l)+br for some integers a,b > 
such that K = a + b is even. Therefore we can represent 
the lattice A as a disjoint union of vertical strips, A = 
A\ U . . . U Ak , such that any strip has width (along the x- 
axis) r or r — 1, and the total number of strips K is even, 
see Fig.[5J The assumption Eq. fTS")) implies that any non- 
trivial logical operator P G C(S)\S is supported on at 
least two strips. Therefore, Lemma [T] implies that we can 
multiply P with some product of generators S a to clean 
out any chosen even strip Ai, A<±, . . . , Ak- However, a 
support of any generator S a has width at most r and 




FIG. 2: Theorem]]]; sketch of the proof. Choose any logical 
operator P G C(S)\S. Assuming di(S) > r, we can clean out 
any vertical strip of width r. Repeating this procedure on 
parallel disjoint strips (even strips), we obtain an equivalent 
logical operator P' which has support on the remaining set 
of disjoint strips (odd strips). If the strips are wide enough, 
then P' commutes with <S on every odd strip individually. It 
must be that the restriction of P' on some odd strip is not 
in iS, otherwise P' G 5. Thus we obtain a non-trivial logical 
operator whose support is contained in a single vertical strip 
which contradicts the assumption that di(S) > r. 

thus it overlaps with at most one even strip, so we can 
apply Lemma [1] repeatedly to clean out all even strips. 
In other words, we can define a new logical operator 

P' = PS, for some S = JJ J] 

p=2,4,...,K aeJ(A p ) 

such that P' is supported only on odd strips, that is 

P' = P[Pi---P' K _ 1 , P;e?(4 3 ). 

Next we note that each generator S a overlaps with at 
most one odd strip. Therefore the inclusion P' G C(S) 
is possible only if Pj G C(S) for all j = 1, 3, . . . , K - 1. 
Since P 1 ^ S there must exist at least one odd strip 
Aj such that Pj ^ S, that is, Pj is a non-trivial logical 
operator. However it contradicts Eq. (|16p since the linear 
dimension d\(Pj) is at most r. □ 

Theorems [1] and [2] are direct corollaries of Proposi- 
tion [TJ Indeed, choose any logical operator P G C(S)\S 
with a linear dimension di(P) < r. Then the weight of 
P is at most rL ^ 1 , that is, d < rL ^ 1 which proves 
Theorem [TJ Furthermore, if D = 2 then the bound 
di(P) < r implies that the support of P can be cov- 
ered by a quasi-lD vertical string. One can construct a 
walk 7 G W(J, P) that implements P in a row-by-row 
fashion. At any intermediate step the contribution to 
the energy cost of a partially implemented P comes only 
from the two end-points of the string. The assumption 
that any qubit participates in 0(1) generators S a implies 
that e m ax{l) = 0(1). It proves Theorem[5J 

Remark: After having tried to construct 3D stabi- 
lizer codes with d = 0(L 2 ) which would saturate the 
bound in Eq. ([2]), we conjecture that any 3D stabilizer 
code with local generators satisfies d = O(L). In fact, 
all known examples of stabilizer codes correspond to a 
bound d = 0(L L D / 2 J ). The pattern that seems to emerge 
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is that there is a trade-off between the minimum weight 
of, say, the logical X versus logical Z operator. The 
ID repetition code shows that one can have a logical 
X of weight O(L), but then the logical Z is of weight 
O(l), so the distance of the code is O(l). In 2D, one can 
make both operators of weight O(L) as the surface codes 
demonstrate. In 3D, the 3D toric code is an example 
where again one logical operator is of weight 0(L 2 ) but 
the other logical operator is of weight 0(L). In 4D, both 
operators are of weight 0(L 2 ). 

III. SUBSYSTEM CODES 

A. Bounds on The Code Distance 

Let us start by generalizing the cleaning lemma to sub- 
system codes. 

Lemma 2 (Cleaning Lemma for Subsystem Codes). Let 
Q be a gauge group of a subsystem code and M be an 
arbitrary subset of qubits. Then one of the following is 
true: 

(1) There exists a non-trivial logical operator P £ C{S)\Q 
whose support is contained in M, 

(2) For any logical operator P G C(Q) one can choose a 
stabilizer S G S such that PS acts trivially on M . 

Proof. Let (il) C V be a subgroup of phase factors. By 
definition of the subgroups Sm and Q(M) one has the 
following inclusion: 

(il)-g(M) QC(S M )nV(M). (17) 

If the inclusion Eq. (jTTJ) is strict, there exists a Pauli 
operator P G V(M) such that P G C{S M ) but P £ 5, 
that is, P G C{S)\Q. It corresponds to case (1) in the 
statement of the lemma. Otherwise Eq. (fTT|) is an equal- 
ity, that is, (il) ■ G(M) = C(S M ) n V(M). Taking the 
centralizer of both sides one gets 

S M =C(g(M))r\V(M). (18) 

Let P G C(g) be any logical operator and Pm be a restric- 
tion of P onto M. By definition, P M G C(G(M))r\V{M). 
It follows from Eq. (|18p that Pm G Sm, that is, one can 
find a stabilizer S G S such that PS acts trivially on 
M. □ 

Note that cleaning out the region M may require 
changing the operator P globally (not only in a small 
neighborhood of M) since the stabilizer group S of a 
subsystem code may lack a local generating set. This is 
the reason why for subsystem codes we cannot clean out 
multiple disconnected regions simultaneously: cleaning 
out one region by multiplying with elements in S may 
generate support in another 'far-away' region which we 
already cleaned out. 

On the other hand, if the stabilizer group S = 5nC(5) 
has local generators, S = (Si, . ■ ■ , S m ), the operator S in 



Lemma [2] can be chosen such that it includes only those 
generators S a whose support overlaps with M. In this 
case the cleaning can be performed "locally" . 

For subsystem codes we shall need a new lemma which 
relates the distance of a code g defined on the entire lat- 
tice A to the distance of a code 5 m obtained by restricting 
g onto some subset of qubits M C A (here we use nota- 
tions introduced in Section [TT|. Recall that Zoo-distance 
between vertices u,v G A is defined as maxj^ \uj —Vj\. 

Lemma 3 (Restriction Lemma). Let g = (Gi, . . . , G m ) 
be a gauge group of a subsystem code on a D -dimensional 
lattice A. Suppose the support of any generator G a can 
be bounded by a hypercube with r D vertices. Let M C A 
be an arbitrary subset and dM be a subset of vertices 
u G A\Af such that loo-distance between u and M is at 
most r. Consider the subsystem code with gauge group 
gM ■ Then one of the following is true: 

(1) The code gM has no logical qubits, 

(2) The code gM has distance at least d — \dM\ 
where d is the distance of g defined in Eq. |P|). 

Remark: By abuse of notation we regard the code 5m 
as a code defined on qubits of M only (otherwise any 
Pauli operator on a qubit u G A\M would be a logical 
operator for the code 5m)- One can always extend the 
code 5m to the rest of the lattice in a trivial way, for 
example, by adding a gauge operator Z u for any qubit u G 
A\M. For simplicity we shall ignore these technicalities. 

Proof. Indeed, suppose 5m has at least one logical qubit 
and let a" be the distance of 5m- Then there exists a 
logical operator P G C(5m) H V(M), P $ 5m and a 
gauge operator G G 5a/ such that \PG\ — d'. Then 
P G C(5)\5, that is, P is a non-trivial logical operator 
for the original code. By definition of the group 5m, 
we can extend the operator G G 5m beyond M to some 
gauge operator G' G 5 such that G is the restriction of G' 
onto M. Since 5 has local generators, such an operator 
G' can be chosen to have support only in M and in dM. 
It means that \PG'\ < d' + \dM\. Since PG' is a non- 
trivial logical operator for the code 5, its weight must be 
at least d, or d < d' + \dM\. □ 

Now we are ready to apply these Lemma's to prove 
Theorem [3] For simplicity we present a proof for D = 1 
and open boundary conditions. Generalization to higher 
dimensions and periodic boundary conditions is straight- 
forward. In ID we have to prove the upper bound 

d < 3r. (19) 

Let M C A be the smallest contiguous block of qubits 
such that the subsystem code 5a/ obtained by restricting 
5 onto M has at least one logical qubit (if there are 
several minimal blocks M, choose anyone of them). Let 
d! be the distance of the code 5ai ■ The restriction lemma 
implies that 

d' >d- 2r. (20) 
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Let us assume that d! > r and show that it leads to a con- 
tradiction. Indeed, d! > r implies \M\ > r, so we can par- 
tition M into three contiguous blocks, M = ABC, such 
that \B\=r-l and A, C ^ 0. Let P, Q e C(Q M )\g M 
be any pair of anti-commuting logical operators for the 
code Qm ■ Since the size of B is smaller than the distance 
of Qm, we can apply Lemma [2] to the code Qm to clean 
out the region B which yields a pair of "cleaned" logical 
operators 

P,Q eC(g M )\G M , PQ = -QP 

such that P, Q act trivially on B, that is, P — PaPc, Q — 
QaQc- Obviously, PQ = -QP iff P A Q A = ~Q a Pa or 
PcQc = —QcPc- Assume without loss of generality 
that P A anti-commutes with Qa (otherwise, relabel the 
blocks A and C) . Since the code Qm has local generators 
of size at most r and the block B has length r — 1, no 
generator of Qm overlaps with both A and C. Therefore 
P, Q G C(Qm) implies Pa,Qa G C(Qm)- Combining this 
observation with the fact that Pa and Qa anti-commute 
we conclude that neither of Pa, Qa can be an element of 
Qm, that is 

Pa, Qa G C(Q m )\Qm: PaQa = -QaPa- 

Consider the subsystem code Qa obtained by restricting 
Qm onto A. We want to argue that Pa, Qa are non-trivial 
logical operators for the code Qa- Indeed, by definition, 
Pa, Qa G C(Q M ) implies Pa, Qa € C(Qa)- Since P A and 
Qa anti-commute, neither of them can be an element of 
Qa, that is, 

Pa,Qa£C{Q a )\Qa- 

Summarizing we have found a contiguous block A C M , 
\A\ < \M\, such that the code Qa has at least one logical 
qubit. This however contradicts the minimality of M and 
thus proves that d! < r. Taking into account Eq. ([20)) we 
arrive at d < 3r. 

The proof of Theorem [3] for 2D subsystem codes is 
analogous. In that case we would consider the minimal 
width strips M of height L such that Qm has at least one 
logical qubit. The same approach generalizes to arbitrary 
dimensions D. 

B. Energy Barrier Upper Bounds 

Let us start by describing some general properties 
of the Hamiltonian H defined in Eq. ([5]). Recall that 
we define a gauge group Q C V as a group gener- 
ated by Pauli operators G a in the Hamiltonian Eq. {SJ. 
We define a stabilizer group S as the center of Q, i.e., 
S = Q DC(Q). The stabilizer group induces a decompo- 
sition of the Hilbert space of n qubits into ^-invariant 
sectors, (C 2 )®" = (J) s £ s , where s labels different syn- 
dromes (irreps of S). By definition, H is block-diagonal 
with respect to this decomposition, 



Each sector C s possesses a subsystem structure C s — 
Ciogicai, s ® £gauge,s- The Hamiltonian H s acts trivially 
on the 2 fc -dimensional logical subsystem Ci og icai,s and has 
a non-trivial spectrum of gauge qubit excitations (which 
may be different for different sectors s). 

The question of whether such Hamiltonian can be a 
model of a self-correcting memory is more involved than 
for stabilizer codes and goes beyond this paper. One of 
the first questions is to determine whether the ground- 
space of the Hamiltonian is confined to one sector, and 
if so, whether the excitations from this sector, call it sq, 
to the other sectors are gapped. The fulfillment of such 
conditions would show that the Hamiltonian has some 
stability with respect to perturbations that affect the log- 
ical qubits. This condition is naturally obeyed for stabi- 
lizer code Hamiltonians. Note that we do not require the 
gauge-qubit excitations in a fixed sector to be gapped, 
since these excitations do not harm the state of the log- 
ical qubits. In particular, there may be low- lying or 
gapless gauge-qubits excitations from the ground-space 
which give rise to a local order parameter, see for example 
the discussion of the quantum compass model in 0] . The 
presence of a local order for the gauge qubits, i.e. a local 
operator acting on the gauge qubits which distinguishes 
two different ground-sectors, does not at all prevent the 
possibility of a 'topological order' for the logical qubits. 
It only expressed the fact that the gauge qubits are not 
a good place to put the quantum information. 

Let us define the gap 

A S = \(H S )-\{H S0 ), (21) 

where X(H S ) is the smallest eigenvalue of H s and Sq is 
the sector that contains the smallest eigenvalue of H. 

A proper definition of the energy barrier that must 
be traversed in order to perform a logical operation is 
dependent on the energy spectrum of the Hamiltonian 
in different s sectors. Consider a logical operator E G 
C{S)\Q and let 7 G W{I,E) a walk on the Pauli group 
implementing E. In order to define the energy cost of a 
Pauli operator E, we consider what the operator does on 
the s sector. Any state in the s sector will be mapped 
onto a state in some other fixed sector s. We define the 
energy cost of E as 

e(E) = A s , £(so)=s . (22) 

Let the maximum energy barrier along a path be 
£max(l) = maxp £7 e(P). Then we define the energy bar- 
rier for E as 

d*(E) = min e max {-f). 

7 6W(I,B) 

The minimum over these energy barriers for any logical 
operation is then given by 

d x = min d x (E). (23) 
Eac(s)\g 



s 



In the definition of d and we minimize over elements 
of C(S)\Q which are comprised of a 'bare' logical operator 
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(an element of C(Q)\Q) times any gauge operator (an 
element of Q), see Eq. (flU)) . The minimum weight of 
operators in C(S)\Q may be much less than the minimum 
weight of 'bare' logical operators. A good example is 
the Heisenberg model on a D-dimcnsional lattice where 
Q is generated by nearest-neighbor XX, YY and ZZ 
operators. If L is odd, the operators X a \\ = Jl MgA X u and 
Z a \i = Il«eA Z u commute with Q but are not generated 
by Q. In addition we have S = I . This implies that the 
minimal weight of the bare logical operators in C(G)\G is 
L D . However, the distance d of this code is 1, because we 
can reduce the weight of, say, X a \\ to 1 by multiplying 
with it with XX-type gauge operators. If some gauge 
qubits operations are costly to implement (i.e. require 
creating high-energy excitations) then it may be more 
physically reasonable to not count these operations for 
free in minimizing the distance and the energy barrier. 
Whether such choices are warranted however depends on 
the features and spectrum of the Hamiltonian and cannot 
be analyzed using the coding framework alone. 

In order to set the stage for proving no-go results which 
do not depend on detailed properties of the Hamiltonian, 
we will derive a simple upper bound on the energy cost 
of a Pauli operator E defined in Eq. (|25J). We apply this 
analysis to the Hamiltonian defined in Eq. © assuming 
for simplicity that r a = — 1, that is, 

m 

H = -Y / G a . (24) 

0=1 

Let E £ V be any Pauli error mapping the ground 
state sector C So to some sector C s and let \ipo) be a 
ground state of H So , such that X(H So ) = (t/'oI-HIV-'o}- 
Since E\tpo) G C s we get an upper bound X(H S ) < 
{i) Q \E^HE\i) Q ). It implies 

e(E) = X(H S ) - X(H S0 ) < (ip \E^HE - H\^ ) 

and thus 

e(E)<\\E^HE-H\\<\\[H,E}\\. 

We can bound the norm of the commutator [H, E] taking 
into account that E either commutes or anti-commutes 
with any generator G a - It yields 

e(E) < 2 #{a : EG a = -G a E}. (25) 

Summarizing we can bound the energy cost of E by 
(twice) the number of terms in the Hamiltonian anti- 
commuting with E. Note that this upper bound is iden- 
tical to the definition of energy cost for stabilizer codes, 
Eq. Q. 

The upper bound Eq. ([23]) opens up the possibility 
for a generalization of the energy barrier arguments in 
Section [H] to subsystem codes. For example, one can use 
Eq. (gU) to show that S =0(1) for the 2D Bacon-Shor 
code. This easily follows from the fact that a partially 
implemented logical operator (say, X — Yii^,i) OIU y 



anti-commutes with the gaug e op erators ZZ at the end- 
points. (Numerical results of [2l| indicate that d* decays 
exponentially with the lattice size L.) 

Using Eq. (j25|) we can also show that d* — 0(d). 
Indeed, choose a non-trivial logical operator P = 
P Ul ...P Ud G C(S)\Q with weight d. The number of 
generators G a anti-commuting with any single-qubit op- 
erator P Ui is at most 0(1). Therefore the number of 
generators anti-commuting with P is at most 0(d). Im- 
plementing the sequence of the single-qubit Pauli oper- 
ators P Ul , . . . , P Ud in an arbitrary order one gets a walk 
7 G W(I, P) with 

1 -max (7) = 0(1). 

C. Proofs of Theorem 1* and 2* 

Let S — (S\,...,S m ) be a stabilizer code with lo- 
cal generators, and let X%, Z\, . . . , X g , Z g be logical 
Pauli operators on some subset of g logical qubits (the 
ones we want to treat as gauge qubits). Let Q = 
(S, X\, Z\, . . . , Xg, Zg). We can regard Q as the gauge 
group of a subsystem code. Clearly, G^C(Q) = S, so S is 
the stabilizer group of Q. In order to prove Theorems 1* 
and 2* let us define a linear distance of a subsystem code. 

Definition 2. Given a subsystem code Q define a linear 
distance of Q as 

di(G) = mm di(P). (26) 
Pec(s)\g 

Recall that di (P) is a linear dimension of P, that is, the 
minimal length of a contiguous interval R C {1, . . . , L} 
such that i?x{l,...,L} D ~ 1 covers the support of P, see 
Section [ill Thus a linear distance of a code character- 
izes the minimal linear dimension of non-trivial logical 
operators along some fixed coordinate axis. 

Proposition 2. Let Q be a subsystem code on a lattice 
A = {1, . . . , L} D with open or periodic boundary condi- 
tions. Assume that a stabilizer group S = Q n C(Q) has 
local generators. S = (Si, . . . , S m ), such that the support 
of any generator S a can be covered by a hypercube with 
r D vertices. Then 

d 1 (Q)<r (27) 

for any L > 2(r — l) 2 . 

Proof. The proof is almost identical to the proof of 
Proposition [T] in Section [ill It goes by assuming that 
d\(G) > r and showing that it leads to a contradiction. 
The only difference is that instead of cleaning out a logi- 
cal operator P £ C(S)\S using Lemma[T]one has to clean 
out a logical operator P £ C(Q)\Q using Lemma Af- 
ter cleaning out all even strips A 2 , . . . , Ak one gets an 
equivalent logical operator P 1 £ C(Q)\Q that has sup- 
port only on odd strips Ai,A 3 , . . . , Ak-i (we used the 
notations from the proof of Proposition [T]) . Since any 
generator of S overlaps with at most one odd strip, the 
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inclusion P' G C(Q) C C(S) implies that a restriction of 
P' onto any odd strip is an element of C(S). At least one 
of these restrictions is not an element of Q. Therefore, 
a restriction of P' onto some odd strip is an element of 
C(S)\Q. This is in contradiction with the assumption 
that di(<?) > r. □ 

Theorems 1* and 2* are direct corollaries of Proposi- 
tion O Indeed, choose any logical operator P € C(S)\Q 
with a linear dimension d\{P) < r. Then the weight of 
P is at most rL 0-1 , that is, d(Q) < rL ^ 1 which proves 
Theorem 1*. Furthermore, if D = 2 then the bound 
di(P) < r implies that the support of P can be covered 
by a quasi-lD vertical string and thus one can construct 
a walk 7 e W(J,P) with e ma x{l) = 0(1). It proves 
Theorem 2*. 



IV. DISCUSSION AND OPEN PROBLEMS 

This paper addressed the problem of constructing 
a self-correcting quantum memory based on stabilizer 
codes with geometrically-local generators. We devel- 
oped several technical tools for analyzing such codes and 
proved upper bounds on the distance that are tight for 
ID and 2D codes. In addition, we defined an energy bar- 
rier separating different logical states and proved that 
for 2D stabilizer codes the energy barrier can not grow 
with the lattice size. We note that a similar conclusion 



has been independently reached by Kay and Colbeck [23| 
using completely different techniques. 

It would be interesting to prove our conjecture that 
d = 0(L) for stabilizer codes in 3D. In addition, we would 
like to bound the energy barrier eP for 3D stabilizer codes 
and for 2D subsystem codes. It might be possible that 
a self-correcting quantum memory can only be based on 
3D subsystem codes such as the one suggested by Ba- 
con [l|. This would require a detailed analysis of the 
energy barrier for these systems. 

Another interesting problem concerns algorithms for 
computing the distance of stabilizer codes with local gen- 
erators. The dynamic programming algorithm allows one 
to compute the distance of stabilizer (subsystem) codes 
with local generators in time of order Lexp(L I?_1 ). It 
requires exponential time for D > 2. It remains an open 
question whether the distance can be computed in poly- 
nomial time for D > 2. 



V. ACKNOWLEDGEMENTS 

We thank Dave Bacon for posing some of the questions 
that lead to the results in this paper. We thank Panos 
Aliferis for many discussions on the topic of this paper. 
We are grateful to Roger Colbeck for explaining to us 
several technical points of [23j . 



[1] D. Bacon, Phys. Rev A 73, 012340 (2005). 

[2] A. Kitaev, Annals Phys. 303, 2 (1997). 

[3] S. D. Sarma, M. Freedman, and C. Nayak, 

Phys. Rev. Lett 94, 166802 (2004). 
[4] A. Kitaev, Annals Phys. 321, 2 (2005). 
[5] L. Jiang, G. K. Brennen, A. V. Gorshkov, K. Hammerer, 

M. Hafezi, E. Demler, M. D. Lukin, and P. Zoller, Nature 

Physics 4, 482 (2008). 
[6] S. Gladchenko, D. Olaya, E. Dupont-Ferrier, B. Doucot, 

L. B. Ioffe, and M. E. Gershenson, Superconducting 

nanocircuits for topologically protected qubits (2008), 

arXiv:0802.2295. 
[7] P. Bonderson, M. Freedman, and C. Nayak, 

Phys. Rev. Lett 101, 010501 (2008). 
[8] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, J. 

Math. Phys. 43, 4452 (2002). 
[9] Z. Nussinov and G. Ortiz, Phys. Rev. B 77, 064302 

(2008). 

[10] C. Castelnovo and C. Chamon, Phys. Rev. B 76, 184442 
(2007). 

[11] R. Alicki, M. Fannes, and M. Horodecki (2008), 

arXiv:0810.4584. 
[12] A. Kay (2008), arXiv:0807.0287. 
[13] D. Gottesman (1997), arXiv:quant-ph/9702029. 
[14] S. Bravyi and A. Kitaev (2008), arXiv:quant- 

ph/9811052. 

[15] H. Bombin and M. A. Martin-Delgado, Phys. Rev. Lett 
97, 180501 (2006). 



[16] H. Bombin and M. Martin-Delgado, Phys. Rev. Lett 98, 
160502 (2007). 

[17] H. Bombin and M. Martin-Delgado, Phys. Rev. B 75, 

075103 (2007). 
[18] E. B. Davies, Comm. Math. Phys. 39, 91 (1974). 
[19] R. Alicki, M. Fannes, and M. Horodecki, J. Phys. A: 

Math. Theor. 40, 6451 (2007). 
[20] D. Poulin, Phys. Rev. Lett 95, 230504 (2005). 
[21] J. Dorier, F. Becca, and F. Mila, Phys. Rev. B 72, 024448 

(2005). 

[22] P. Aliferis and A. Cross, Phys. Rev. Lett 98, 220502 
(2007). 

[23] A. Kay and R. Colbeck (2008), arXiv:0810.3557. 

[24] Note that the Ising models are in fact not stable against 
perturbations (stray magnetic fields etc.) while Kitaev's 
toric code model is stable against sufficiently weak per- 
turbations. 

[25] The assumption that logical qubits are encoded into the 
ground-state of the system may fail for some models such 
as topological quantum computation where information 
is encoded into a collective state of anyonic excitations 
rather than the ground-state. 

[26] Our analysis can be easily generalized to Hamiltonians 
H — —^2™ =1 r a Sa with arbitrary real coefficients r a - 

[27] Such dephasing in the energy eigenbasis can be explicitly 
enforced by measuring the stabilizer generators S a using 
an ancilla qubit. We do not need to learn the outcome of 
this measurement, since the role of this measurement is 



12 



only to project onto a certain eigenvalue of the stabilizer 
[28] A generalized toric code on a D-dimensional lattice has 
qubits living on hypercubes of dimension D = \ L D/2\,X- 
type stabilizers living on hypercubes of dimension D' + 1, 
and Z-type stabilizers living on hypercubes of dimension 
D' — 1. Logical X-type operators correspond to closed hy- 
persurfaces of dimension D' and logical Z-type operators 
correspond to closed hypersurfaces of dimension D — D' 
on the dual lattice. 
[29] If we would only consider the minimum weight of ele- 
ments of C(Q)\Q, that is the weight of 'bare' logical op- 
erators on the logical qubits, one can easily construct ID 
counterexamples to Theorem 1*. Take (odd) n copies of 



a small code, such as the Steane code [[7,1,3]] and define 
two high-weight logical operators as X = X1X2 ■ ■ -X n 
and Z = Z1Z2 ■ ■ ■ Z„. Here Xi, Zi are the logical oper- 
ators for the ith Steane code. It is clear that even when 
we multiply these operators X and Z with elements in S, 
their minimum weight will scale with n. Aside from this 
logical qubit there are also gauge qubits. If we multiply 
X, Z with the logical operators of these gauge qubits, we 
can reduce the weight of X or Z to O(l). 
[30] It is well known that the codes corresponding to different 
choices of the syndrome s are equivalent up to a local 
change of basis in every qubit. 



